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Abstract 

We evaluate binomial series with harmonic number coefficients, providing 
recursion relations, integral representations, and several examples. The results 
are of interest to analytic number theory, the analysis of algorithms, and cal- 
culations of theoretical physics, as well as other applications. 
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Introduction 

The evaluation of harmonic number sums has been useful in analytic number the- 
ory for some time (e.g., [21 E]). Recently, the evaluation of Euler sums [51 [10] has 
been important in various areas of theoretical physics, including in support of Feyn- 
man diagram calculations (e.g., [7]). More recently, the performance of generalized 
harmonic number sums has been useful in evaluating Feynman diagram contributions 
of perturbative quantum field theory [8]. In addition, harmonic number sums often 
arise in the analysis of algorithms. This especially applies to algorithms involving 
searching, sorting, or permutations (e.g., [HI [18]). 

Here we are interested to evaluate sums containing simultaneously two types of 
special numbers of enumerative combinatorics: binomial coefficients and harmonic 
numbers. Our methods complement those of Refs. [3, [T5|. The analytic approach 
of [5] is based upon hypergeometric summation and relies on identities of generalized 
hypergeometric functions p+i-Fp at the very special argument of 1. These identities 
include the Chu-Vandermonde formula for 2-^1, the Pfaff-Saalschiitz theorem for 3F2, 
the Dougall-Dixon theorem for 5F4, and the Whipple transformation for yFg. In 
Ref. |i5j, symbolic computation using the Newton- Andrews-Zeilberger algorithm 
was used to discover harmonic number relations. These relations were motivated by 
considerations to prove certain 'super congruences' for Apery numbers. 

Our results not only complement those of Ref. [3 [T5|, but demonstrate that vari- 
ous generalizations are possible, including the introduction of a summation parameter. 
The latter feature means that many previous results may be considered special cases, 
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and points out that many more results should be achievable in the future. 

After first introducing some notation, we show how various recursion relations for 
the sums of interest may be developed directly. We then illustrate the development 
of integral representations for these sums, and provide examples. To emphasize the 
usefulness of the integral representations, we present low order cases explicitly. In 
addition, we give a representation of binomial-harmonic number sums in terms of the 
generalized hypergeometric function and its derivative. In particular, with regard to 
the Gauss hypergeometric function, we are then able to develop results in terms of 
Legendre function Pj,. 

In the final section, we introduce a variation. Namely, we consider binomial har- 
monic sums over the order of the generalized harmonic numbers. Additional con- 
siderations in analytic number theory motivate the study of these and related sums 
[9]. In particular, a certain combination of these sums represents a truncation of a 
dominating sum for the Li/Keiper constants. 

We put Hn = Y,k=i the usual harmonic numbers, and H^''^ for the gener- 

alized harmonic numbers 



^r^E^:' Hr^^H^'l (0.1) 



These are given in terms of polygamma functions ip^^'^ as 

= [^^^~'\n + 1) - ^^^-'Hl)] , (0.2) 

where ^/'^'""-'^^(l) = (— l)'"(r — l)!C(r) and ( is the Riemann zeta function. The asymp- 
totic form of Hn is well known, if„ = Inn + 7 + o(l), where 7 is the Euler constant. 



Indeed, by Euler-Maclaurin summation we have 



2n Jn 



°° Pi(x) , 
—dx, 



n 



(0.3) 



where the periodic Bernoulh polynomial Pi (a;) = Bi{x — [x]) = x — [x] — 1/2. The 
asymptotic form of Hj{^ for large n is immediately known from that of ip^^'~^''{n + 1), 
and we have (e.g., p. 260) 



^(")(^) = (-1) 



n-l 



(n-l)\ n\ ^/ 1 
^ + t: — rr + 



n+2 



(0.4) 



We investigate here a subclass of generalized harmonic number sums 



\z\ < 1. 



(0.5) 



In particular, in this paper we restrict for the most part to m = 1. When, in addition, 
g = r = 1, we write 



n ' 



\z\ < 1. 



j=0 \J/ 
When simply p = 0, we omit the superscript. 

We note a very convenient recursion for the general sums of Eq. (0.5): 

d 

S^f^^^ (g, r, m, z) = ^-q^S^^ ^ ^) ■ 



(0.6) 



(0.7) 



In this way, from an initial sum, successive sums may be obtained. 

Although we do not follow this line of inquiry, we mention that many integral 
representations for binomial coefficients are known. These (e.g., [Tl], pp. 372-375; 
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p. 53) include [II], (p. 375) 



2n+2 ^7r/2 2"^^ f^l'^ 

cos" X sin nx sin 2mx = / cos"" x cos nx cos 2mx dx. 



i n ' 

\mj TC Jo TT Jo 

(0.8) 

As a contour integral, we have for complex a and integral j > 1, 

— / (1 + 2)"2-^-M2, 0<r<l. (0.9) 



\j J 2Tii 

We may note that this equation is particularly well formed for taking derivatives with 
respect to the parameter a. If both n > 1 and j > 1 are integral, we no longer have 
a branch point at z = —1, and may write 



<r <oo. (0.10) 



Equations (0.9) and (0.10) may be immediately verified by using the binomial theorem 
to compute the residue of the integrand at z = 0. 

Recursion relations for binomial-harmonic number sums 

We have 

Proposition 1. We have the recursion relation (a) 

5„+i(l) = 2^„(l) + ^-— , (1.1) 

n + 1 



and (b) 



(z + l)"+i — 1 

5„+i(^) = (1 + z)S„(z) + ^-—^^ . (1.2) 

n + 1 



For part (a), 5*1(1) = 1, and for part (b), Si{z) = z. 



Proposition 2. We have the recursion relation (a) 



^Si(l) = 25(^)(1) + E (^^j SiP + 2 E ^+i^^(2, 2, . . . , 2, 1 - 1, 1, . . . , 1, 3; -1), 

(1.3) 

and (b) 



p-1 



i=0 



.iF,(2,2,...,2,l-n;l,l,...,l,3;-z) 

(1.4) 



For part (a), 5'i^''(l) = 1, and for part (b), Si'\z) = z. 
Proposition 3. We have the recursion relation (a) 



S^^\z)=nzSn-i{z) + {l + zr-l. 



(1.5) 



Let 



P4p, z) = J2 f (■_,)- = ^ p-iFp-2i2, 2, . . . , 2, 1 - n; 1, 1, . . . , 1; -z). (1.6) 
j=i \3 ^/ J 



Then we have (b) 



5^)(z) = .ze7^;^^-(^) 



5;ii(z) +n(3n{p,z). 



(1.7) 



Proof of Proposition 1. The proof of part (b) subsumes that of part (a). By using 
a recursion relation for the binomial coefficient [IjJ (p. 822), we have 



n+l 



Sn+l{z) = 



+ 1 



j=o 



n+l 
j=0 



n\ n 



SJz, 



(1.8) 



where we have used the facts = — Hq. Further, by shifting the summation 

index and using the recursion relation for harmonic numbers, we have 



1 



j=0 \ I ^ ^ 



n ' 



J + 



1 



n + 1 



[1.9) 



wherein the latter sum may be obtained by integrating the binomial theorem. 
Proof oi Proposition 2a. For this part, omitting the z — 1 argument, we have 



n+l 

i=o 



n+l 
3=0 



J. 



+ 



' n ' 
,J - h 



n+l 



3=1 



' n 

KJ - K 



j=0 



[1.10) 



We next interchange the order of the two sums, and to complete this part of the 
Proposition we need to perform the sum 



E 



J 



n ' 



{ ij + 1) \J. 



E 

j=0 



(2), 
(1).- 



n 



j + 2\j + ir 



(1.11) 



where {a)j — r{a+j) /T{a) is the Pochhammer symbol, and P is the Gamma function. 
We further use (2)^/(3)^ = 2/(j + 2) and 



n 



;i.i2) 



Appealing to the series definition of e+iFe completes part (a). The steps for part (b) 
are very similar, and are omitted. 

Proof oi Proposition 3a. We use the property — n(^z]): so that 



'n-l\ 1 



+ Hi-, 



n 



n 



,=1 V^' - V \J 



n— 1 



n — 1 



n — l\ 



3 

n-2 



n — 1 



n — 1\ z^ 



■1=1 



J 



n 



n-l 

E 



;;>'-"'-""---"^""--"^|:(;:J)7 



n 



zSn_,{z) + -{{i+zr-i) 



n 



(1.13) 



wherein the latter sum may be obtained by integrating the binomial theorem. 
For part (b), we proceed similarly, 



j=l \J J 
^nf:f-'(''~^]H,z^ 



n 



n-l 



Er._ljH,z^+n^-'Hr.z- 



= n 



p-1 



i=i 



+ H.i z^ + nP-^H„ 



n 



n—l 
j=Q 



(1.14) 
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We then binomially expand, interchange sums, and the Proposition is complete. 

Remark. Similar considerations can be given for quadratic and more generally 
nonlinear sums, that are outside of the present investigation. 

Examples 

The recursion relations (1.1) and (1.2) with initial condition may be explicitly 
solved. For Eq. (1.1) we have 



5,(1) = 2"if„ + ^— 2Fi(l,l;n + 2;-l)-2"ln2. 



n + 1 



:i.i5) 



Since easily we have 



n + 



we obtain 



1 "1 
— 2Fi(l, l,n + 2- -1) - 2^^ In 2 = -2" J2 J^-. 



^„(l) = 2"|i7„-f:-i-,. 



:i.l6) 



:i.l7) 



a known result [15]. For Eq. (1.2) we find 

1 1 



SJz) = il+zT 



H„ 



(;z + l)"+i (n+1) 



2F1 (l,r2 + l;n + 2;^^^ +ln^ ^ 



z + lj \z+l 

(1.18 



By using a transformation formula [TT] (p. 1043), we may write this as 

S4z) = (1 + zr + , ]^., 2F1 (1, 1; n + 2; --) + In 

z[z + lp[n + l) \ zj \z+l 



:i.i9) 



For z = 1 here we have the explicit reduction to Eqs. (1.15) and (1.17). At z = —1 
we have 

(1.20) 



.^aM;^ + 2;i)^I^±fM.!i±l. 

1 "^[n + 1) n 
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Therefore, we obtain 1) = —1/n. 
From Eqs. (1.5) and (1-18) we obtain 



S^^\z) = nz{l + z)""-^ 



^n-i + 7— -tt:::- ai^i f 1, ri; n + 1; 



(2; + 



z + 1 



In 



z+1 



(1.21) 



Remark. Prom the recursion relations of Proposition 2, it appears that the portions 
of the respective binomial-harmonic series containing harmonic numbers are given for 
p > 1 by 



Slf\l) = 2"-^(n), 



n-p 



(1.22) 



and 

S^\z) = (l+^)"-^(n). 



z{z + l)''-P (n+l-p) 



2Fi(l,l;n-p + 2;-^) -In 

(1.23) 



Integral representations for linear binomial-harmonic number sums 

We have 

Proposition 4. We have the integral representation (a) 

S'<^\z) = nz f [t pF^_i(2, 2, . . . , 2, 1 - n; 1, 1, . . . , 1; -zt) 

J 



dt 



pFp_i(2, 2, . . . , 2, 1 - n; 1, 1, . . . , 1; -z)] j—^ 



(2.1) 



and (b) 



(q, 1, 1, ^) = P^P-i(2, 2, . . . , 2, 1 - n; 1, 1, . . . , 1; -zt) 

- pFp_i(2, 2, . . . , 2, 1 - n; 1, 1, . . . , 1; -^)] ^-f—^dt. 



(2.2) 



10 



Proof. For part (a), we employ the relation Hj = iplj + 1) + 7, where ip = F'/F is 
the digamma function, together with an integral representation for this function [TT] 
(p. 943): 



j=0 



The integral is absolutely convergent and we may interchange summation and inte- 
gration. We then rewrite the summation, using relation (1.12), and apply the series 
definition of pFp_i. Equation (2.1) follows. 
For part (b), we use Eq. (0.2), writing 

^^)(g, 1, Iz) = tfl^^'-'^J + 1) - ^^'-'Kl)] (^) 

Here, we employed the result of multiply differentiating an integral representation for 
the digamma function to obtain that for the polygamma function. The integral is 
again absolutely convergent and we may interchange it with the summation. Carrying 
out the summation, we have Eq. (2.2). 

Remarks. The integral representation of the polygamma function used above, 

^(«-i)(^) = ( ^1^) In'^-i t dt = (-l)^(g - l)!C(g, z), (2.5) 



lo \ t 

shows the connection with the Hurwitz zeta function ({s, a) at integer values of q. 
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Reference [8] contains a number of other integral representations for harmonic 
numbers. Similarly, these may be applied to obtain other integral representations of 
the sums SIp\ 

The direct verification of the general property (0.7) from the integral representa- 
tions requires the result 

^pFp_i(2,2,...,2,l-n;l,l,...,l;-zt) 

= (n - 1)2P-H pFp_i(3, 3, . . . , 3, 2 - n; 2, 2, . . . , 2; -zt), (2.6) 

as well as a transformation formula for p-Fp_i. It is possible that this transformation 
may be effected by noting 

(2-n),= + (1-n),. 

(2.7) 

To show the utility of integral representations for binomial-harmonic number sums, 
we specialize in the following two sections to low order cases of the sums S!;^\ and 
provide corresponding details. 

Explicit expressions for sums S^\z) 

We have 

Proposition 5. We have for \z\ < 1, 

Sn{z) = nz{l + zy-^ ,F2 (l, 1, 1 - n; 2, 2; ^) , (3.1) 

giving 

12 



i3l 
(2), 



" ' jV 1 

1 + - = 7 

2 J 2P-1 



(2), 
.(1). 



+ 1 



Corollary 1. 



5«(z) = (i+zr^|i+ 



1 + z 



(3.2a) 



and 



2n-l + 



1 — n 



[l + z)"" 



+ n^z{l + nz) 3F2 ( 1, 1, 1 - n; 2, 2; 

(3.26) 



1 + z 



The right side of Eq. (3.1) is easily verified to be a polynomial of degree n in 
z, as it must. For the (terminating) 3F2 function is a polynomial of degree n — 1 in 
z/{l + z). When multiplied by the prefactor of z{l + z)"'~^, we obtain a polynomial 
of degree n. We could continue Corollary 1 indefinitely, but these instances serve our 
current purpose. 

Proof. As with Eq. (2.3) at p = 0, we have 



Sn{z)-Er,]^' 

j=0 \J J 



n\ , fH^ -1 



Jo t - 1 



dt 



^- f\{l + z)--{l + zt)-]-^. 
Jo 1—1 

Changing variable with u — 1 — t, we have 



(3.3) 











zu 


n ^ 




[ l + z\ 





du 
u 



(3.4) 



and then obtain the Proposition. 

Corollary 1 follows by applying property (0.7). 
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Explicit expressions for some sums ^^^^(l) in terms of harmonic numbers 



We have 
Proposition 6. We have 



Sn{l) = 2" 

^«(l)=n2"-i 



n 1 



»t— 1 -j^ 



1 + 2^ 



(4.1a) 
(4.16) 



and 



in-2 



n(n + l"; 



n-2 ^ 
^n-2 - E - 



+ [2n(l - 2-") - 1] 



2n 
n — 1 



(4.1c) 



Ahhough Eqs. (4.1a) and (4.1b) are known resuhs [15], our method of proof is 
different. In addition, the latter appears to be a rediscovery of earlier closed form 
summations 

Proof. Using the integral representation of the previous section, writing Eqs. (3.3) 
and (3.4) at ^ = 1, we have 



SJl] 



1 

[-2"- (1 + t)"] 



u 



t- 1 
du 



u 



1 (i-wY 
ipin + 1) + 7 + / ^ —dw - In 2 



2" 



'1/2 W 

2Fi(l,l,n + 2;-l) -ln2 



(4.2) 



2^^(71+1) 

Here we have applied Eqs. (A5) and (A6) of Ref. [9J. By Eq. (1.16), Eq. (4.1a) 
follows. 
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Writing Eq. (2.3) a,t z — 1 and p = 1, we have 



^ ^ / -dt 



/■I df 
= n / [t(l + ir-^-2"-M-^ 
7o t — 1 



n-l 



t+ 1 



n-l 



1' 



(4.3) 



Changing variable with u — 1 — we have 

5«(l)=n2-^ 



1-^ 
2 



n-l 



u Jo \ 2 



n— 1 



= n2 



n-l 



n-l 



- E 



2'^(2-" - 1), 



(4.4) 



giving Eq. (4.1b). 

Following a similar procedure, writing Eq. (2.3) at 2; = 1 and p = 2, we have 



n\ , /-i P - 1 



dt 



= n2-2 / 
^0 

Now with u — 1 — t, we have 



j=o \JJ Jo t-1 

/■I df 
= n / [nt(l + nt)(l + t)""' - 2"-2n(n + 1)1— - 
io t — 1 

t+ 1^ n-2 



nt(l + nt) 



n{n + 1) 



(it 



1 



(4.5) 



S(^\l)^2--'\n{l-u)[(n + l)-nu]{l-'l 



n-2 



n{n+l) } — 
u 



2"-" <^ n(n + 1) 



1 - 



u 



n-2 



- 1 



u 
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2"-2 { n{n + 1) 
2 



n-2 

4n 



2 4-77 1 

+ 2n) -(2^-" - 1) + -[1 - 2-"(l + n)] . (4.6) 

n — 1 n — 1 ) 

Simplifying the latter terms on the right side of this equation gives Eq. (4.1c), and 
completes the Proposition. 

Hypergeometric approach 

We have 

Proposition 7. For integers p > 1, 



X 



X d 



pFp_i[-i/, . . . , -I/; 1, . . . , 1; {-lYx\\^=ri 



Proof. We have 



p du 

+HnX-'' pFp_i[-n, . . . , -n; 1, . . . , 1; {-Ifx] 



so that 



Using 



and the sum 



J2[ .] x^^ . . . , -n; 1, . . . , 1; (-l)^'x]. 



d_ 
dv 



{-uYj = p(-i^)M^ + 1) - ^(^ - 3 + 1)], 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



n-j 



j=Q \3/ 

we obtain the Proposition. 



t(llYHrnX--^X-S(^Ul,p,l,l), (5.5) 
m=0 V'V ^ 
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We focus on the p — 2 case of Proposition 7, and obtain a number of series 
representations. For this, we introduce the Legendre functions of the first kind Pi,, 
and the function 



Rn{z) 



dP.{z) 



du 



_ln('i±iVn(;^), 



(5.6) 



where n > is an integer, and Pn is a Legendre polynomiaL We have 
Proposition 8. 



1,2,1, 



X 



+Hr^x 2Fi{-n,-n;l;x). 
^(1 - x^K (^) + H^x-^{1 - xTPr, 



l + x 



X 



(5.7) 



Proof. We first observe that through the use of a transformation formula for 2-F1 
we have 

P,(l - 2x) = 2Fii-u, i/ + 1; 1; x) = (1 - x f 2F1 (^-u, -u; 1; , (5.8) 



givmg 



n(.)^(if)\..(-.-..i^). 



z-V 



(5.9) 



Differentiation of this equation with respect to u and comparison with the defining 
relation (5.6) yields 



fz + iy d f z-i 

Rniz) = 7— 2F1 —v, —y\ 1: 

With a change of variable, we obtain Eq. (5.7). 



(5.10) 
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We now discuss and illustrate Proposition 8, in the course of which we obtain 
Corollary 2 

2 



5^(1, 2, 1, 1) = J:H, Q = (2if„ - H,^) {^^^ 



(5.11) 

recovering a known result |15j . 

There are several series representations of the function i?„ available, including the 
Bromwich forms [1] 



A 1 



k=l 



n-1 



2k + 1 



H,.(,-)^2g(-ir- („_,„„,,,^ 



[Pfc(^)-P„(z)], 



(5.12) 



(5.13) 



and from a formula of Jolliffe |12j . 



Rjz) = -2\n {^^]Pr,{z) 



1 



[z-lY{z+lY\xi 



z + l 



(5.14) 



Besides these older results, i?„ and closely related polynomials have been very recently 
restudied [17]. Since the Legendre polynomial has PniX) = 1 and Pk{ — 1) = (—1)^, 
we readily see that the function satisfies -Rn(l) = and -R„(— 1) = 2(— l)"iJ„. It 
is evident from the Bromwich formula (5.13) that 



R4z) = 2{H,^ - H^)P4z) + 2 ^(-1)"+'=- TT7^-T—T-Pk{^) 



(5.15) 



The property i?„(— 1) = 2(— l)"iJ„ follows immediately from this equation. With 
the use of the duplication formula satisfied by the digamma function, the property 
i?n(l) = may also be deduced from Eq. (5.15). 



Equivalent to the Chu-Vandermonde summation 2-Fi(— —n; 1; 1) = (^^), is the 
relation, via Eq. (5.9), 

„^^(,_,),.p„(i±£).p. (5,6) 

In addition, with the help of Eq. (5.15) and the use of partial fractions, we have 

\im(l-xrRn(l^)=-2(^'']{H2n-Hn). (5.17) 
2:^1 VI — X/ \^/ 

Therefore, from Proposition 8 we obtain Corollary 2. 

By the same method of this section, on taking more derivatives with respect to u 
in Eq. (5.1), it is possible to get higher values of q in the sums defined in Eq. (0.5). 
In this regard, we very briefly mention the functions 



Qix,z)= ,F,ix,-x;l;z) = J2^-f^z\ Q(x) ^ Qix, -1), {b.li 

j=0 

and the expansion 



Q{x) = 1 + J2A2jX^^, (5.19) 
written in Ref. [16]. Obviously the Maclaurin coefficients here are given by 



(5.20) 



x=0 

Then by the property (5.4) at p = 1, the coefficients A2j involve sums and products 
of generalized harmonic numbers. Further connections with Legendre functions and 
Stirling numbers of the first kind we do not pursue here. 

Among the generalized hypergeometric functions that reduce to polynomials is 
the family 

Hn{x,a,z)= 3F2(-n,n + l,a;; 1, a; z), (5.21) 
19 



where n > is an integer and x, a e C, with a ^ —n — 1, —n — 2, . . . These functions 
may be obtained by a certain integration over shifted Legendre polynomials, and 
similarly for higher order functions pFp_i. Obviously the set of functions of Eq. (5.21) 
includes i?„(n +1,1,2;) and related functions. One may then ask for a generalization 
of Proposition 8 for the case of p = 3 in Proposition 7, an effort that we are leaving 
to future work. 

We also mention an integral representation for the sum of Proposition 8, and how 
the case of Corollary 2 may be otherwise obtained. We have 



5^(1, 2, 1,^)= [\^F,{-n,-n,l,zt)- ^Fii-n, -n,l, z)]-^. 

Jo t — 1 



(5.22) 



As special case, we have 
-5^(1,2,1,1) = 



1 



2Fi{-n, -n, l,t)-l 



'2n^ 
,nj JO 

/2n^ 
\ n , 



1 



(l-a;)-"P,(l-2a;)-l 



Pn{z) - 1 



dt 
t-l 

dx 
x-l 

dz 



l + z' 



(5.23a) 



(5.236) 



(5.23c) 



where we changed variable and used Eqs. (5.8) and (5.9). One way to demonstrate 
equality with the result of Corollary 2 is to show that the integral of this equation 
satisfies the same recursion relation Qn+i — Qn = ^/ (2n + 1) + l/(2n + 2) — 2/(n + 1) 
as the quantity Qn — i?2n — 2if„, with initial condition Qi — —1/2. We also note the 
integral representation 

.1 (^n _ 1)2 



Qn 



/ 

Jo 



-dt. 



(5.24) 
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Other class of binomial harmonic sums 



Finally, we describe another class of sums, where the summation is now over the 
order of the generalized harmonic numbers. 



m=2 



m, 



M 



Z . 



(6.1) 



We have 

Proposition 9. Let be the associated Laguerre polynomial of degree n (e.g., 
[HEIIII]). We have 



1 r 



Sn{M,z)=z Ll_^{z\xit) - n 



t^' -I 
t-1 



dt. 



(6.2) 



Proof. We use relation (0.2) together with an integral representation for the 
poly gamma function. 



5„(M,z)= 



m=2 



\m—l 



m - 1 ! 



^(™-i)(M + l) -V'^"'"^^!) 



E 

m=2 



n 



(_1)-1 rift 



M 



t dt. 



m-l)\Jo \ t-1 



(6.3) 



Interchanging the finite summation with the integration and applying the power series 
definition of Ll^_i, we obtain the Proposition. 

Remarks. In the limit as M — > oo in Eq. (6.2), we have the representation 

dt 



lim 5„(M,z) 

M^oo 



Ll^_i{z\nt) — n 



t-l 



(6.4) 



This is the same limit in which i^l™"* — > C{m). In particular, we put 



Siin)= lim [SniM,-l] 

M^oo 



5„(M,-l/2)] 



(6.5) 
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With a change of variable, we obtain agreement with the integral representation of 
Si{n) of Ref. [9j (p. 216). In particular, this alternating binomial sum provides the 
apparently dominant contribution to the Li/Keiper constants A„ of the Li criterion 
for the Riemann hypothesis [9]. This sum has been shown to be O(nlnn), describing 
the significant cancellation within the summand. 
Proposition 10. We have 



j=2 \ J / 

Proof. By using the recursion relation of generalized harmonic numbers, we have 
from the definition (6.1) 



This recursion may also be found from the integral representation of Proposition 9. 
We form the difference 




(6.6) 




(6.7) 




(6.8) 



We then integrate by parts. 





zn 



M 




(6.9) 
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In the last step, we may put t = exp(— f), and use an extension of the Laplace 
transform of the Laguerre polynomial (e.g., formula 7.414.2 of [TT], with A = —z and 
/x = 0). 

Since we have the initial value ^^(l,^;) = — 1 — n^; + (1 + z)", solution of the 
recursion (6.7) gives the Proposition. 

An alternative proof of Proposition 10 is now obvious. For from binomial expan- 
sion we obtain 



r 



1=0 



J / j=i e=o 

= E ^'4? + M + uzHm. (6.10) 

It is possible to expand the factor {t — 1)"^ of Eqs. (6.2) and (6.4) as a geometric 
series, thereby giving another series representation of the sum Sn-, and providing a 
third proof of Proposition 9. We have 
Corollary 3. We have 



Sn{M,z) 



'UzHm - M + 

j=l L 



1 + 



(6.11) 



A simple shift of summation index here recovers the relation (6.6). 
In order to obtain the result (6.11), we write 



S^{M,z) = -zY: / [Ll_,{-zv)-n][e- 



{]\I+j+l)v _ g-(j+l)t) 



]dv. 



(6.12) 



j=0 



We then use for Re > 0, 



I 

Ll_^i-zv)e-''''dv = - 



(6.13) 
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Remarks. Corollary 3 must be equivalent to using the partial fractions form of 
coming from that for the polygamma function, in (6.1), and then interchanging 
summations. 

By manipulating Sn{M, z) we may obtain related sums. For instance, by integrat- 
ing Eq. (6.1), using Proposition 8, and interchanging integrations, we have 



E 

m=2 



, m / m + 1 



up,{n + l) 



1 + 



u 



n+l 

n+l 



t 



M 



t - 1 



dt 



3 



n 



- M u. 

2 



(6.14) 



From Proposition 10 we obtain the combination 



M 



5„(M,-l)-5„(M,-l/2)=5] 



2J, 



n 



H 



M 



M 

E 

j=l L 



(6.15) 



The latter expression is the M < oo form of Eq. (19) of Ref. [9]. 

Plainly, the sums of this section may be generalized to those such as 



Sn{M,p,q,r,z) = J2 



n 



m=2 



(6.16) 
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